This project displays an interactive land use/land cover and watershed distribution map of the Hawaiian Islands. This map can be utilized to determine the composition of the land use and land cover within a specific watershed. This can be useful for hydrological modeling projects within Hawaiian watersheds.
For Task 3, make a map of land use / land cover and watersheds for the big island of Hawaii.
Familiarize yourself with the data attributes (see link), and download the shapefile data for each. After reading the shapefile data into R, make finalized map(s) in which you:
There is flexibility! You can decide:
Whatever you decide is great, just make sure that your final outputs are presented in a nice professional HTML that you’d be proud to share with someone so that they can see your awesome spatial code!
library(tidyverse)
library(janitor)
library(kableExtra)
library(naniar)
library(skimr)
library(sf)
library(tmap)
library(paletteer)
library(lubridate)
library(here)
watershed <- read_sf(dsn = here(".", "Watersheds"),
layer = "Watersheds") %>%
st_transform(crs = 6628) # Set the coordinate reference system
lulc <- read_sf(dsn = here(".", "Land_Use_Land_Cover_LULC"),
layer = "Land_Use_Land_Cover_LULC") %>%
st_transform(crs = 6628) %>%
filter(landcover != 0) # Remove land use/land cover category of 0
hawaii_clip_tmap <- tm_basemap("Esri.WorldImagery") + # Put in a basemap
tm_shape(watershed) +
tm_borders("lightblue", lwd = .5) +
tm_shape(lulc) +
tm_layout(legend.outside = TRUE, legend.outside.position = "right") +
tm_fill("landcover", title = "Land Cover", palette = c("orange", "purple", "blue", "yellow"), alpha = 0.5) +
tm_layout("Land Use/Land Cover for Watersheds of the Hawaiian Islands")
tmap_mode("view") # Sets it to interactive view
hawaii_clip_tmap